Field-induced barrier transparency of Bloch waves in tight-binding lattices 
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A rectangular potential barrier for a Bloch particle in a tight-binding lattice is shown to become 
fully transparent by the application of a strong ac field with appropriate amplitude and frequency. 
Such a curious phenomenon bears some connection with the field-induced barrier transparency 
effect known for freely-moving particles scattered by an ac-driven rectangular barrier; however, 
for a Bloch particle transparency is not related to a resonant tunnneling process across the cycle- 
averaged oscillating potential barrier, as for the freely-moving quantum particle. The phenomenon 
of field-induced transparency is specifically discussed here for photonic transport in waveguide arrays 
and demonstrated by full numerical simulations of the paraxial (Schrodinger) wave equation beyond 
the tight-binding approximation. 
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I. INTRODUCTION 



Coherent control of electronic, photonic or matter wave 
transport in driven semiconductor superlattices, quan- 
tum dots, waveguide arrays and optical lattices has re- 
ceived a great and continuous interest over the past two 
decades [TH3], stimulating a wide number of experimen- 
tal and theoretical investigations in different physical sys- 
tems [SUlOj. Examples of coherent control by the appli- 
cation of ac fields include, among others, the coherent de- 
struction of tunneling between two wells in a bistable po- 
tential [5] and the suppression of quantum diffusion and 
dynamic localization in tight-binding lattices [5]. Such 
phenomena have been demonstrated in a series of recent 
experiments [5141 Oj. mainly based on transport of light 
waves or cold atoms in waveguide arrays or optical lat- 
tices, where dephasing and many-body effects can be ne- 
glected. Other interesting phenomena occur in the tun- 
neling and scattering processes of ac-driven free particles 
across potential wells or barriers [TT1I14) . In particular, a 
freely-moving quantum particle can be resonantly trans- 
mitted across a periodically-driven rectangular potential 
barrier, in spite the probability of tunneling through the 
static potential barrier is almost zero HT, ^V2} . This phe- 
nomenon, originally predicted for quantum tunneling and 
referred to as field- induced barrier transparency I12j , 
is rather generic and can be similarly observed for opti- 
cal tunneling at modulated dielectric interfaces [15]. In 
the accelerated (Kramers-Henneberger) reference frame, 
the ac-driven barrier behaves like a periodically-shaken 
potential barrier [TTJ [T3j. At high oscillation frequen- 
cies, the oscillating rectangular barrier can be replaced 
at leading order by its cycle-averaged (static) barrier, 
whose profile shows a characteristic double-hump shape 
that sustains metastable states. Thus resonant transmis- 
sion of the incident particles across the ac-driven barrier 
observed at certain below-barrier energies can be simply 
explained as a resonant tunneling process, similar to e.g. 
the Ramsauer-Townsend effect, the resonant electronic 
tunneling across a double-barrier structure, or the reso- 
nant light transmission of light in a Fabry-Perot optical 



cavity. In this work we consider the coherent motion of 
a Bloch particle in a tight-binding lattice with a strong 
impenetrable rectangular potential barrier, and demon- 
strate that the application of an ac field with appropriate 
amplitude and frequency can make the potential barrier 
fully transparent. This phenomenon, which can be again 
referred to as field-induced barrier transparency for Bloch 
particles in analogy to its counterpart for free particles, 
has however a very different origin. In particular, barrier 
transparency is observed for incident Bloch waves at any 
allowed energy in the band, i.e. it is not related to reso- 
nant tunneling of the cycled-averaged potential as in the 
free particle case. Here we investigate the phenomenon of 
field-induced barrier transparency of Bloch waves by con- 
sidering photonic transport in periodically-curved waveg- 
uide lattices [H [U] , however the present analysis is rather 
generic to coherent transport in driven tight-binding lat- 
tices and could be therefore of interest to other physical 
set ups, such as to cold atoms or Bose-Einstein conden- 
sates in accelerated optical lattices JOJ. 
The paper is organized as follows. In Sec. II, the ba- 
sic model is described with specific reference to light 
transport in photonic lattices, and a theoretical analysis 
of field-induced barrier transparency in ac-driven tight- 
binding lattices is presented. In Sec. Ill the theoretical 
predictions obtained in the high-frequency regime are 
confirmed by direct numerical simulations, based on both 
the time-periodic coupled-mode equations of the driven 
tight-binding lattice and the Schrodinger wave equation 
with a one-dimensional potential. Deviations from the 
predictions based on the averaged lattice model in the 
high frequency modulation limit, such as a shift of the 
resonant curves, are briefly discussed. Finally, in Sec. IV 
the main conclusions are outlined. 
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II. FIELD-INDUCED 
BARRIER-TRANSPARENCY IN AC-DRIVEN 
TIGHT-BINDING LATTICES: THEORETICAL 
ANALYSIS 

The starting point of our analysis is provided by a 
standard model of light transport in a one-dimensional 
array of tunneling-coupled optical waveguides with a 
periodically-curved optical axis [11 [9] . In the waveguide 
reference frame, light transport is described by the fol- 
lowing Schrodinger-type wave equation for the electric 
field envelope (t){x,z) (see, for instance, [H I16j) 

(1) 

where A = X/{2tt) is the reduced wavelength of light 
waves, rig the substrate refractive index, x and z are the 
transverse and the longitudinal spatial coordinates, re- 
spectively, and V{x) is the optical potential, which is re- 
lated to the refractive index profile of the straight array 
by the simple relation V{x) ~ ns—n{x). The last term on 
the right hand side of Eq.(l) is a fictitious refractive index 
gradient arising from axis bending and with a z-varying 
slope given by F{z) = —nsXo{z), where xo{z) is the axis 
bending profile and the dot indicates the derivatives with 
respect to 2: [1]. In its present form, the paraxial wave 
equation (1) is formally equivalent to the Schrodinger 
equation describing the dynamics of a quantum parti- 
cle of mass rig in the potential V{x) driven by a time- 
dependent external force F(z), provided that the spatial 
z variable in the optical system is replaced by the tem- 
poral variable in the quantum problem, and the photon 
wavelength A is replaced by the Panck's constant. The 
potential V is given by the superposition of the periodic 
and barrier potentials, i.e. V{x) = Vp{x) + Vb{x), where 
Vp{x + a) = Vp{x) is the periodic potential describing 
a homogeneous array of equally-spaced wells, and Vb{x) 
is the barrier potential, which is assumed to describe a 
rectangular barrier of height Vq and width L ^ a [see 
Fig. 1(a)]. The role of the barrier potential is to intro- 
duce a step, of height Vq, into a finite number of wells in 
the lattice, say from the waveguide n = to the waveg- 
uide n = A^, as shown in Fig. 1(a). The potential height 
Vq is typically assumed to be much larger than the width 
of the lowest energy band of the array, yet smaller than 
the gap between the first and second lattice bands, as 
shown as an example in Fig. 1(c). In this way, in the 
absence of the ac driving force any Bloch wave packet, 
belonging to the lowest band of the lattice and traveling 
along the lattice, is reflected from the potential barrier 
[see Fig. 1(c)]. To study the role of the ac driving force 
and the possibility to make the barrier transparent, let us 
introduce the nearest-neighboring tight-binding approx- 
imation and let us assume that the lattice is mostly ex- 
cited in its lowest-order Bloch band. Such conditions are 
satisfied, for example, by considering an array of weakly- 
coupled waveguides which is initially excited by a broad 
beam tilted at an angle smaller than the Bragg angle 
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FIG. 1. (a) Behavior of the optical potential V{x) — ns—n{x) 
(solid curve) of a one-dimensional lattice with a barrier, com- 
posed by the superposition of the periodic optical potential 
(lattice period a) and of the rectangular barrier potential 
(height Vo, width L, represented by the dashed curve in the 
figure). The dotted curve in the figure is the intensity profile 
of the incident Gaussian wave packet in the numerical simula- 
tions of Sec.III.B. (b) Band diagram (three lower-order bands) 
of the periodic optical potential. The energy in the vertical 
axis is defined as the eigenvalue of the Hamiltonian Ho of 
the periodic potential entering in Eq.(l). (c) Space-energy 
band diagram of the lattice with the rectangular barrier. The 
height Vo of the barrier is larger than the bandwidth 4A of 
the lowest lattice band, but smaller than the gap separating 
the first two bands, (d) Behavior of the Bloch excitation co- 
efficients versus incidence angle 9 oi & plane wave, normalized 
to the Bragg angle 63 = A/ (2a). 



Sb = A/(2a) (HUrnj. Under such assumptions, from Eq. 
(1) the following coupled-mode equations can be derived 

mm 

iCn = -A(c„+i + C„_i) - f{z)nCn + (JPnCn (2) 

for the amplitudes c„ of the field trapped in the individ- 
ual waveguides, where A > is the coupling constant 
between adjacent waveguides. 



f{z)^lF{z) = -"^i,{z) 



(3) 
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is the normalized forcing, 



where we have set 



and 



1 for < n < iV 
otherwise . 



(4) 



(5) 



To study field-induced barrier transparency, it is worth 
introducing, in place of c„, the amplitudes a„ defined by 
the relations 



in / dUiCi + if^PnZ 



OLn = c„ exp 

so that Eqs.(2) take the form 

m„ = -A„(z)a„+i - A*_i(z)a„_i 
where we have set 



An{z) = A exp 







(6) 



(7) 



(8) 



Let us now assume a sinusoidal bending of waveguide axis 
with spatial frequency oj = 27r/A and amplitude A, i.e. 



so that 



where 



xo{z) — y4cos(27rz/A), 
f{z) = /oCos(27rz/A), 



A2A 



(9) 
(10) 

(11) 



Let us also assume that: 

(i) the barrier height Vq and the quanta of modulation Xuj 
are much larger than the width 4AA of the tight-binding 
energy band; 

(ii) the modulation frequency is chosen such that the res- 
onant condition 



lu) = (T 



X 



(12) 



is satisfied for some (small) integer / (typically I — 1 
OT I = 2). In the high-modulation frequency limit [as- 
sumption (i)], at first order in a multiple-scale asymp- 
totic analysis of Eqs.(7) the spatial evolution of the am- 
plitudes a„ is dominated by the cycle-average coupling 
rates A„ {1/ A) dzAn{z) (see, for instance, [l7]). 
i.e. one has 



+ 1 



A* 



(13) 



Taking into account that exp[zr sin(a;2;)] = 
J„(r) exp(ma;z) and using Eqs.(5) and (10), 
from Eq.(8) it readily follows that 

AJo(r) for n ^ -l,iV 
A„ = <^ AJ,(r) for n = -l (14) 
(-l)'AJ/(r) for n = N 



r ; 



/o _ 2-KansA 
uj XA 



(15) 



Equations (13) thus describe light transport in an effec- 
tive lattice with two defects in the coupling rates at lat- 
tice sites n = —1 and n = N. However, if the amplitude 
A of modulation is tuned such that 



Jo(r) = ±J,(r), 



(16) 



the effective lattice is homogeneous (i.e. defect-free), and 
thus any Bloch wave packet propagates in the lattice 
without being reflected [IS]. That is, if the modulation 
frequency and amplitude are tuned to satisfy the condi- 
tions (12) and (16), the effect of the ac driving is to make 
the barrier Vb fully transparent for Bloch wave packets 
belonging to the lowest energy band. 



III. NUMERICAL SIMULATIONS 

In this section we confirm, by direct numerical simu- 
lations of both the tight-binding equations (2) and the 
full wave equation (1), the field-induced barrier trans- 
parency phenomenon predicted in the previous section 
in the high-frequency modulation regime. 



A. Tight-binding lattice model 

Let us first consider beam reflection and transmis- 
sion in the tight-binding lattice model described by the 
coupled-mode equations (2) with periodic coefficients. 
Equations (2) have been numerically integrated using an 
accurate fourth-order variable-step Runge-Kutta meth- 
ods with absorbing boundary conditions. As an exam- 
ple, Figs.2(a-d) (upper panels) show the numerically- 
computed evolution of lattice site occupation probabil- 
ities |c„(2;)p for a/A — 4, N = 6, lu = a and for a few 
increasing values of the modulation amplitude, measured 
by the dimensionless parameter F = /o/w. For compari- 
son, the lower panels in Figs.2(a-d) show the correspond- 
ing evolution of occupation probabilities as obtained by 
integration of the averaged effective lattice equations 
[Eqs.(13) and (14)], valid in the high modulation fre- 
quency limit. The array is initially excited by a Gaus- 
sian distribution c„(0) — exp[— (n — ng)'^ / w^] exp{iqn) 
with mean position no — —26, width w = 4 and mo- 
mentum q = Tr/2 (corresponding to the largest cycle- 
averaged group velocity far from the barrier). Note that 
at F = 1.435, at which the condition (16) is satisfied for 
1 = 1, the wave packet is not reflected, and the barrier ap- 
pears to be fully transparent [see Fig. 2(c)]. In the simula- 
tions shown in Fig. 2, the ratio w/A = 4 is large enough to 
ensure the validity of the averaged equations (13), at least 
at first-order approximation. In fact, a more careful com- 
parison of the results obtained from the original coupled- 
mode equations with periodic coefficients [Eqs.(2) or (7)] 
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FIG. 2. Beam evolution (snapshots of lattice site intensities |c„(2)|^) for an initial Gaussian wave packet distribution with 
transverse momentum q as obtained by numerical simulations of the coupled-mode equations (2) (upper panels) and by the 
averaged equations (17) (lower panels) for a/A — 4 and oj = a. In (a-d), the initial transverse beam momentum is g = 7r/2, in 
(e) q = 7r/3, whereas in (f) q = 7r/4. In (a) F = (non-modulated lattice), in (b) F = 0.8, in (c),(e),(f) F = 1.435 (transparency 
condition), in (d) F = 1.8. 
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FIG. 3. Behavior of the power reflection (i?) and transmission 
(T) coefficients of a Gaussian wave packet at the first barrier 
discontinuity near the first transparency resonance bj = a 
numerically computed (a) using the tight- binding model (2), 
and (b) the full wave equation (1). Parameter values in the 
two cases are given in the text. 



and the averaged equations (13) shows some slight dis- 
crepancies, which basically arise from neglecting higher- 
order terms in the asymptotic analysis of Eqs.(7) (for 
more details see, for instance, [ITj). In particular, ac- 
cording to the average model (13), transparency should 
be observed regardless of the initial value of wave packet 
momentum g, i.e. of its mean energy; however, numer- 
ical simulations of the original (periodic) coupled-mode 
equations (2) show that some reflected light is observed 
at the barrier when the initial wave packet momentum q 
(i.e., its mean group group velocity) is reduced. This 
is shown, as an example, in Figs. 2(e) and (f), where 
the evolution of lattice site occupation probabilities are 
depicted for the same initial Gaussian wave packet of 



Fig. 2(c) (i.e. at the transparency condition), but with 
an initial momentum lowered to g = 7r/3 [Fig. 2(e)] and 
q = 7r/4 [Fig. 2(f)]. Another difference between the aver- 
age and periodic coupled-mode equations can be seen by 
computing the first resonance curve of the transparency 
process, depicted in Fig. 3(a). The figure shows the be- 
havior of the power reflection R and transmission T coef- 
ficients versus uj (near the first resonance w ~ ct = Vq/A) 
for the same initial Gaussian wave packet with momen- 
tum q = 7r/2, as obtained by numerical integration of 
Eqs.(2). For the sake of simplicity, the coefficients R and 
T have been computed by considering the first poten- 
tial step solely of the barrier, so that multiple reflections 
that would arise in the presence of the two potential dis- 
continuities of the rectangular barrier [see, for example. 
Fig. 2(b)] are avoided . For each modulation frequency a;, 
the modulation amplitude A was correspondingly varied 
such that the ratio F = /o/w [see Eq.(15)] remains con- 
stant and equal to 1.435, at which the condition (16) is 
satisfied. Note that, the condition of exact transparency 
(i? = 0) is attained at the ratio uj/cr ~ 0.95, which is 
slightly smaller than 1, as expected from the averaged 
model [Eq.(12) with I — \\. As discussed above, the va- 
lidity of the averaged model becomes more accurate, and 
the agreement between the averaged and original (peri- 
odic) coupled mode equations closer, as the ratio ~ cr/A 
of barrier height and width of tight-binding energy band 
is increased. 



B. Full-wave equation 

We checked the validity of the tight-binding lattice 
analysis and the onset of field-induced barrier trans- 
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FIG. 4. Beam intensity evolution (snapshot of \(j}{x, z)\'^) of 
an initial Gaussian beam along the periodically-curved waveg- 
uide array of Fig. 1 (a) for increasing values of the modulation 
amplitude [measured by the parameter F, given by Eq.(15)]: 
(a) F = (straight array); (b) F = 0.8, (c) F = 1.435 (trans- 
parency condition), and (d) F = 1.8. In all the simulations, 
the input Gaussian beam is tilted at the angle 6 = 6b/'2., and 
the modulation frequency is a; = 0.98a". 
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FIG. 5. Same as Fig. 4, but for an initial Gaussian beam with 
F = 1.435, Lo = 0.98a and tilt angle (a) 6 = 63/2,, and (b) 
e = 6»s/4. 



parency in photonic waveguide arrays by direct numer- 
ical simulations of the paraxial wave equation (1) using 
standard pseudospectral methods. The optical potential 
V{x) — Us — n{x) of the lattice used in numerical simula- 
tions is shown in Fig. 1(a) and corresponds to a typical 
effective index profile of Lithium-Niobate waveguide ar- 
rays, fabricated by the proton-exchange technique and 
probed at A = 1.44 fim wavelength (see, for instance, the 
experiment reported in Ref.|16j). The refractive index 
profile n^{x) of each waveguide in the array is taken to 
be given by [T5] 



nwix) = An 



erf [(a; -t- Wg)/D.j;] - erf [(a; - Wg)/Dx] 
2eri{wg/D,) 



(17) 



where Wg = 3.5 /im is the channel width, = 1 /^m the 
diffusion length, and An = 0.003 the peak index change. 
The waveguide spacing (lattice period) is a = 12 /im, 
and the substrate refractive index at the probing wave- 
length is Us = 2.1381. The band structure of the lattice, 
numerically computed by a standard plane-wave expan- 
sion method, is shown in Fig. 1(b). From the width of 
the lowest order (tight-binding) lattice band, a coupling 
rate A ~ 3.25 cm~^ can be estimated between adja- 
cent waveguides. When the lattice is excited by a broad 
wave packet, tilted at an angle 9, it generally breaks up 
into the superposition of different wave packets, belong- 
ing to the various lattice bands and refracting at differ- 
ent angles (see, for instance, [161 llll)- The fractional 
excitations of the different lattice bands are given by 
the Bloch-wave excitation coefficients C„, defined as in 
Ref.|19j. Figure 1(d) shows the numerically-computed 
behavior of C„ {9) for the various lattice bands of the pe- 
riodic part of the optical potential of Fig. 1(a) versus the 
tilt angle 6 of the incident beam, in units of the Bragg 
angle 9b — A/(2a). As one can see, for broad input 
beams tilted at an angle 9 smaller than half of the Bragg 
angle, the lowest-order lattice band is mainly excited, 
which indicates that in this case the tight-binding model 
of Sec. II can be safely applied. The barrier height used 
in the numerical simulations is Vq = An/5 = 6 x 10"'', 
whereas its width is L = 72 fj,m. A schematic of the 
space-dependent band diagram of the lattice, shown in 
Fig. 1(c), clearly indicates that, in the absence of the 
external ac field, the barrier is impenetrable and any 
Bloch wave packet, incident onto the barrier, is re- 
flected. This is shown, as an example, in Fig. 4(a), where 
the numerically-computed evolution of the field inten- 
sity \(l){x, z)p along the non-modulated lattice is depicted 
by assuming, as an initial condition, the tilted Gaussian 
beam 0(a:, 0) — eyip[—{x + xo)^/w^]exp{2TTi9x/X) with 
spot size w = 30 /im, offset xq = 108 /im [see the dot- 
ted curve in Fig.l(a)] and tilt angle 9 = 9b/2~ 1.719°. 
Let us now introduce a sinusoidal modulation of the axis 
bending xq{z) at a spatial frequency w = 2ti/K close 
to the first resonance [/ = 1 in Eq.(12)], which corre- 
sponds to a spatial modulation period A = 2.45 mm. 
Note that the ratio between modulation frequency lo and 
coupling rate A of waveguides turns out to be ~ 7.9, i.e. 
the high frequency modulation condition is well satisfied. 
Figures 4(b-d) show the evolution of the same Gaussian 
wave packet as in Fig. 4(a) but in the periodically-curved 
waveguide array for increasing values of the modulation 
amplitude A, measured by the dimensionless parameter 
r given by Eq.(15). Note that the cycle-averaged trans- 
verse group velocity of the wave packet, far from the bar- 
rier region and at the tilting angle 9 — 9 3/2, is equal to 
Vg = 2Aa|Jo(r)|, and thus decreases as F is increased 
from zero [Fig.4(a)] to F = 1.8 [Fig.4(d)]. As the modu- 
lation is increased, the Gaussian wave packet is less and 
less reflected from the barrier, till a nearly reflectionless 
regime, corresponding to full barrier transparency, is at- 
tained when the transparency condition (16) is reached 
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[see Fig. 4(c)]. Barrier transparency is observed at dif- 
ferent values of the initial tilt angle 9 (i.e. initial wave 
packet momentum), expect for small tilt angles at which 
the wave packet shows a more complex dynamics at the 
barrier crossing (see Fig. 5). This behavior, already no- 
ticed in Sec. II in the framework of the tight-binding anal- 
ysis, is mainly ascribable to a discrepancy between the 
average model [Eq.(13)] and the original coupled-mode 
equations with periodic coefficients [Eqs.(2) or (7)]. A 
slight discrepancy can be also seen when computing the 
resonance curve of the field-induced barrier transparency 
process using the full wave equation (1), similarly to 
what already noticed within the tight-binding model [see 
Fig. 3(a)]. In Fig. 3(b), the numerically-computed power 
reflection (i?) and transmission (T) coefficients versus uj 
(near the first resonance u) a — Vb/A) are depicted 
for the Gaussian wave packet of Fig. 4. Like for the res- 
onance curve computed using the tight-binding model 
[Fig. 3(a)], the coefficients R and T have been calculated 
by considering the first potential step solely of the bar- 
rier, and for each modulation frequency w the modula- 
tion amplitude A was correspondingly varied such that 
r = /o/w — 1.435. Note that, as in Fig. 3(a), the reso- 
nance frequency turns out to be slightly smaller than the 
theoretical value a predicted by the averaged model. 
As a concluding remark, it is important to stress that, 
as compared to the phenomenon of field-induced bar- 
rier transparency for a freely-moving quantum particle 
[llj in which particle transmission occurs solely at spe- 
cial values of initial momentum (energy) that match the 
metastable states of the cycle-average potential barrier, 
the phenomenon of barrier transparency for Bloch wave 
packets predicted in this work is relatively insensitive to 
the initial wave packet momentum and can not thus be 
ascribed to a resonant tunneling phenomenon as for the 
free particle. However, as opposed to the case of Ref.[TT], 
where the frequency and amplitude of the ac field may 
take relatively arbitrary values |20| . in our system the 
modulation frequency as well as the modulation ampli- 
tude should satisfy certain resonance conditions, namely 



Eqs.(12) and (16). 



IV. CONCLUSIONS 

In this work it has been theoretically shown that a 
rectangular potential barrier for a Bloch particle in a 
tight-binding lattice can be made fully transparent by the 
application of a strong ac field with appropriate ampli- 
tude and frequency. As this phenomenon bears some con- 
nection with the field-induced barrier transparency phe- 
nomenon previously predicted for freely-moving quan- 
tum particles scattered by an ac-driven potential barrier 
pn IT?| . the transparency effect for the Bloch particle 
has a rather different physical origin, as discussed in this 
work. In particular, in the high-frequency limit trans- 
parency is attained independently of the energy of the 
Bloch wave packet, and therefore particle transmission 
across the barrier can not be explained as a resonant 
tunneling process across the cycle-averaged potential, as 
for a freely moving particle . This phenomenon could 
be experimentally observed in periodically-curved waveg- 
uide arrays, where light transport along the lattice mim- 
ics the coherent temporal evolution of a Bloch particle 
with an external ac driving field [9 . Numerical simula- 
tions based on the paraxial (Schrodinger) wave equation 
have also shown that the phenomenon of field-induced 
barrier transparency persists beyond the tight-binding 
approximation, provided that the lattice is initially ex- 
cited by a broad beam tilted at an angle smaller than the 
Bragg angle. 
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